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Abstract 

We discuss certain analogies between quantization and discretization of clas- 
sical systems on manifolds. In particular, we will apply the quantum dy- 
namical entropy of Alicki and Fannes to numerically study the footprints of 
chaos in discretized versions of hyperbolic maps on the torus. 
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1. Introduction. 



Classical chaos is associated with motion on a compact phase-space with high 
sensitivity to initial conditions: trajectories diverge exponentially fast and never- 
theless remain confined to bounded regions [TJEl 



In discrete times, such a behaviour is characterized by a positive Lyapounov 
exponent log A, A > 1, and by a consequent spreading of initial errors 5 such 
that, after n time-steps, 5 t— > 5 n ~ 5 \ n . Exponential amplification on a compact 
phase-space cannot grow indefinitely, therefore the Lyapounov exponent can only 
be obtained as: 

log A := lim — lim log ( -p- 

t->oo n s^o \ o 
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that is by first letting <5 — > and only afterwards n — > oo. 

In quantum mechanics non-commutativity entails absence of continuous tra- 
jectories or, semi-classically, an intrinsic coarse-graining of phase-space determined 
by Planck's constant h: this forbids 5 (the minimal error possible) to go to zero. 
Thus, if chaotic behaviour is identified with log A > 0, then it is quantally sup- 
pressed, unless, performing the classical limit first, we let room for 6 — > |3]. 

In discrete classical systems, one deals with discretized versions of continuous 
classical systems or with cellular automata 0EIIZ| with finite number of states. 
In this case, roughly speaking, the minimal distance between two states or con- 
figurations is strictly larger than zero; therefore, the reason way log A is trivially 
zero is very much similar to the one encountered in the field of quantum chaos, 
its origin being now not in non-commutativity but in the lack of a continuous 
structure. Alternative methods have thus to be developed in order to deal with 
the granularity of phase-space (BHI303IH1III • 

An entropic approach is likely to offer a promising perspective. For suffi- 
ciently smooth classical continuous systems, the exponential spreading of errors is 
equivalent to a net entropy production, better known as Kolmogorov dynamical, 
or metric, entropy (Sj. The phase-space is partitioned into cells by means of which 
any trajectory is encoded into a sequence of symbols. As times goes on, the rich- 
ness in different symbolic trajectories reflects the irregularity of the motion and is 
associated with strictly positive dynamical entropy [TO] . 

A frequency approach to the numerical evaluation of the entropy production 
has recently been applied in discretized version of various chaotic continuous clas- 
sical dynamical systems [TT]. In this paper, we suggest a different strategy: mo- 
tivated by the similarities between quantization and discretization of continuous 
classical dynamical systems, we propose to use the quantum dynamical entropy 
recently introduced by Alicki and Fannes [HUH], which we shall refer to as ALF 
entropy. 

The ALF-entropy is based on the algebraic properties of dynamical systems, 
that is on the fact that, independently on whether they are commutative or not, 
they are describable by suitable algebras of observables, their time evolution by 
linear maps on these algebras and their states by expectations over them. 

As such, the ALF-entropy applies equally well to classical and quantum sys- 
tems, and reduces to the Kolmogorov entropy in the former case. In particular, it 
has been showed that it allows a quite straightforward calculation of the Lyapounov 
exponents of Arnold cat maps on D-dimensional tori [14J. 

In this paper we aim at showing how the ALF-entropy may be of use in a 
discrete classical context, precisely when 2-dimensional cat maps are forced to live 
on a square lattice with spacing -k (N integer), with particular focus upon the 
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emergence of the continuous behaviour when N i — > oo. 

In quantum mechanics, the classical limit is achieved when h — ► 0. Analo- 
gously, in the case of discrete classical systems, by letting the minimal distance 
between states go to zero, one might hope to recover a well defined continuous 
dynamical system, perhaps a chaotic one. Also, very much as in the semi-classical 
approximation, one expects the possibility of mimicking the behaviour of discrete 
systems by means of that of their continuous limits and vice versa. However, this 
is possible only up to a time tb, called breaking-time jlj; it can be heuristically 
estimated as the time when the minimal error permitted, 5, becomes of the order 

of the phase-space bound A. Therefore, when, in the continuum, a Lyapounov 

1 A 

exponent log A > is present, the breaking-time scales as tb = -log — . 

log A o 

In the following, we shall consider discrete dynamical systems obtained by 
discretizing a subclass of the Unitary Modular Group of (2-dimensional) Toral 
Automorphisms |3j containing the well known Arnold Cat Maps. We shall provide: 

• the algebraic setting for the continuous limit N i — > oo; 

• the technical framework to construct the ALF-entropy and numerical short- 
cuts to compute it; 

and study: 

• the behaviour of the entropy production and how the breaking-time r# is 
reached in hyperbolic systems; 

• the differences in behavior between hyperbolic and elliptic systems; 

• the distribution of eigenvalues of the multitime correlation matrix used in 
computing the ALF-entropy; 

• the behaviour of the entropy production in the case of Sawtooth Maps [H3 
[TU[T2j which are discontinuous on the 2-dimensional torus. 



2. Automorphisms on the Torus 

Usually, continuous classical motion is described by means of a measure space X, 
the phase-space, endowed with the Borel cx-algebra and a normalized measure /i, 

n(X) = 1. The "volumes" n(E) = / dfi(x) of measurable subsets E C X represent 

Je 

the probabilities that a phase-point x G X belong to them. By specifying the 
statistical properties of the system, the measure \x defines a "state" of it. 
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In such a scheme, a reversible discrete time dynamics amounts to an invertible 
measurable map T : X \— > X such that \ioT = /i and to its iterates {T^}j & %. Phase- 
trajectories passing through x G X at time are then sequences {T^x}j e z [31 • 

Classical dynamical systems are thus conveniently described by triplets (X, /i, T); 
in the following, we shall concentrate on triplets (X,fj,,T a ), where 

X = T 2 = R 2 /Z 2 = {x = (xi, x 2 ) (mod 1)} (la) 

r «&)=C« a ;)(::) <™ di >- ™ 

dfi(x) = dxi dx2 ■ (lc) 
Remarks 2.1 



i. Since det (T a ) = 1, the Lebesgue measure defined in (flcj) is invariant 
for all a G Z; 

ii. The eigenvalues of ( x + a \ ) are (a + 2 ± ^/(a + 2) 2 -4) /2. They are 
conjugate complex numbers if a G [—4,0], while one eigenvalue A is 
greater than 1 if a ^ [—4,0]. In this case, distances are stretched 
along the direction of the eigenvector |e + ), S a \e + ) = A|e + ), con- 
tracted along that of |e_), S a \e-) = A _1 |e_). 

For such a's all periodic points are hyperbolic [T7| . 

iii. Tt = (fi) is the Arnold Cat Map 0. Then, T x G {T a } aeZ C 
SL 2 (T 2 ) C GL 2 (T 2 ) C ML 2 (T 2 ) where ML 2 (T 2 ) is the subset of 
2x2 matrices with integer entries, GL 2 (T 2 ) the subset of invertible 
matrices and SL 2 (T 2 ) the subset of matrices with determinant one. 

iv. The dynamics generated by T a G SL 2 (T 2 ) is called Unitary Modular 
Group [5 (UMG for short). 



For future comparison with quantum dynamical systems, we adopt an algebraic 
point of view and argue in terms of classical observables, precisely in terms of 
complex continuous functions / on X = T 2 . 



These functions form a C* algebra Ax = C° (X) with respect to the topology 



given by the uniform norm 



o = SU P 



fix) 



the Lebesgue measure \i defines a state on Ax which evaluates mean 
values of observables via integration: 



/ ^(/) := / dx f( x ) 
J x 



(2) 
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• the discrete-time dynamics T a : X i— > X generates the discrete group of 
automorphisms 9^ : Ax >— ► Ax, given by 

e£ (/)(») = (*)), jez, (3) 

that preserve the state, ^00^ = u^. 
Definition 2.1 

The dynamical systems (X,fj,,T a ) will be identified by the algebraic 
triplets (Ax, Up, 6 Q ). 



2.1. "Weyl" discretization 

In the following we shall proceed to a discretization of the systems introduced in 
the previous Section and to the study of how chaos emerges when the continuous 
limit is being reached. 

Roughly speaking, given an integer N, we shall force the continuous classical 
systems (Ax, w^, Q a ) to live on a lattice L N C T 2 given by: 

p E (Z/iVZ) 2 } , (4) 

where (Z/iVZ) denotes the residual class (mod N). 

A good indicator of chaos in continuous dynamical systems is the metric en- 
tropy of Kolmogorov [3j (see Section El below). We can compare discretization of 
classical continuous systems with quantization; in this way, we can profitably use a 
quantum extension of the metric entropy which will be presented in Section HJ To 
this aim, we define a discretization procedure resembling Weyl quantization [T%fTn| ; 
in practice, we will construct a *morphism j7)v,oo from Ax = C° (X) into the abelian 
algebra D N 2 (C) of N 2 x iV 2 matrices which are diagonal with respect to a chosen 
orthonormal basis {\£)}i & (%/N' Il f- The basis vectors will be labeled by the points 
of a square grid of lattice spacing with ^ ^ ^ (N — 1) (N identified with 0) 
superimposed onto X = T 2 . 

In order to define Jn, oo, we use Fourier analysis and restrict ourselves to the 
*subalgebra Wexp £ Ax generated by the exponential functions 

W(n)(x) = exp(27ri n ■ x) , (5) 

where n = (n\,n-2) G Z 2 and n ■ x = n\ X\ + ri2 xi- The generic element of W exp 
is: 

f{x) = Lw(n)(x) (6) 

nel? 




Quantum Dynamical Entropies in Discrete Classical Chaos 7 



with finitely many coefficients f n = JJ ^ x f( x ) e 2mnx different from zero. 
On W exp , formula ([2} defines a state such that 

u,(W(n)) = 5 nfi - (7) 

Further, since n • (T a x) = (T* r n) • x, the automorphisms Q map exponentials 
into themselves: 

Q a (W(n)) = W(T^n), = ( 1 + a f\ (8) 



Remark 2.2 



The latter property no longer holds when a^Zas will be the case in 
Section lOl where we deal with Sawtooth Maps [T5llT6| fT7j . 

Following Weyl quantization, we get elements of D N -i out of elements of Wexp by 
replacing, in (JHj), exponentials with diagonal matrices: 

W(n)^W(n):= ^ e*%* \£) (£\ , l=(^h)- (9) 

£e(Z/7VZ) 2 



Definition 2.2 

We will denote by the *morphism from the *algebra W exp into 

the diagonal matrix algebra D N 2 (C), given by: 

W exp 3 f .— > ^(J) := £ f n W{n) 

= E < 10 > 

Remarks 2.3 

i) The completion of the subalgebra Woo with respect to the uniform 

f(x) is the C* algebra A x = C° {X) [201 - 



norm ||/|| = sup 



ii) The *morphism J§ : W exp ^ ^ (C) is bounded by || 



N.oo 



Using the Bounded Limit Theorem [20], J'ffoo can be uniquely ex- 
tended to a bounded linear transformation (with the same bound) 
Jn,oo '■ -Ax 1 — > D N 2 (C). 

iii) Jn,oo(Ax) =D N 2(<C). 
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We go back from D N z (C) to a *algebra of functions on X by defining a *morphism 
Joo,n that "inverts" Jn,oo in the N — > oo limit. In the Weyl quantization the "in- 
verting" *morphism is constructed by means of coherent states \/3 («)), x 6 X, 
with good localization properties in X. 

Definition 2.3 

We will denote by Joo,n '■ D N 2 (C) i — > Ax the *morphism defined by: 
D N 2 (C) 3 M i — > Joo,n(M) (x) := (f3 (x) \ M \ (3 (x)} , (11) 
where \ (3 (x)) are coherent vectors in TC N 2 = <D n2 

We now construct a suitable family of \(3 (x))\ we shall denote by |_-J and (■) the 
integer and fractional part of a real number so that we can express each T 2 as 
x = (i^i^pf) + ^fe^i^iy Then we associate x E T 2 with vectors of 
Hn 2 as follows: 



x 



|/3(aj)) = An (x)\[N Xl \,[Nx 2 \) + 

+ A 12 (x) | [Nx x \ , LiVa; 2 J + 1) + A 21 (x) \ [Nx x \ + 1, LiVa; 2 J ) + 

+ A 22 (x)\[Nx 1 \ +l,[Nx 2 \ +1) 



(12) 



We choose the coefficients A,-, so that 



x 



1 and that the map (jl2|) be 



(13) 



invertible: 

'A n (x) = cos (f (N Xl )) cos (f (Nx 2 )) 
^ A 12 (x) = cos (f (Nx 1 )) sin (f (Nx 2 )) 
* A 2 i (x) = sin (f (JVari)) cos (§ (JVar 2 )) 

^A 22 (aj) = sin (f (Nxi)) sin (f (Nx 2 )) 

Therefore, from definitions 12.21 and 12.31 it follows that, when mapping Ax onto 
D N 2 (C) and the latter back into Ax, we get: 



f N (x) := (Joo >N o J N ) (f) (x) 



'J(x)\£) 



cos (nfi (Nxi)) cos (tcu (Nx 2 )) 

[N Xl \ +p [Nx 2 \ + a 



x/ 



N 



N 



(14) 



Remarks 2.4 
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i) From ([Till . / = /at on the lattice points. Moreover, although the first 
derivative of (jHj) is not defined on the latter, its limit exists there 
and it is zero; thus, we can extend by continuity to a function in 
C 1 (T 2 ) that we will denote again as /n- 

ii) We note that Ran( l 7 o,Ar) is a subalgebra strictly contained in Ax; 
this is not surprising and comes as consequence of Weyl quantization, 
where this phenomenon is quite typical |I81I19|. 



We show below that J7oo,tv ° J~n,oo approaches \j^ x (the identity function in Ax) 
when N — > 00. Indeed, a request upon any sensible quantization procedure is to 
recover the classical description in the limit K — ► 0; in a similar way, our discretiza- 
tion should recover the continuous system in the — > limit. 



Theorem 1 Given / G A* = C° (T 2 ) , lim 

N— >oo 



{Joo,N i7iV,oo — J-A*) (/) 



Proof: Since A" = T 2 is compact, / is uniformly continuous on it. Further, 



\Nx\ (Nx) , 
denoting cc = — H — — , ^ (jvaj) < 1 implies 

therefore, for all e > there exists Nf >e such that 



[iVa:J 



N 



0; 



iV > jV>, e 



/(*)-/ 



N 



£ 

< - 

2 



uniformly in x. Moreover, according to Remark 12.41 (L). /jy G C°(T 2 ), thus the 



previous inequality holds for /jv, too. Since /at = / on 
sufficiently large iV, 

[iV:z;J 



/ (*) - In (x) 
uniformly in x. 



f(x)-f 



N 



+ 



jV 



N 

[Nx\ 
N 



it follows that, for 



^ e 



3. Kolmogorov metric entropy 

For continuous classical systems (X,/j,,T a ) such as those introduced in Section [2J 
the construction of the dynamical entropy of Kolmogorov is based on subdividing 
X into measurable disjoint subsets {Ei} £=1 2 D such that {J e Et = X which form 
finite partitions (coarse grainings) S. 

Under the dynamical maps T a in (fTbl) . any given £ evolves into T~i(£) with 
atoms T~ J (Ei) = {x G X : T 3 a x G Ei}; one can then form finer partitions £[o,n-i] 
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whose atoms E^^...^ := E io f]T a l {E il ) ■■■f > \T a n+1 (E in _ t ) have volumes 

AW<^i := V (E l0 n T-\E H ) • • • n T-" +1 (^_j) . (15) 
Definition 3.1 

We shall set i — {io i\ ■ ■ ■ i n -i\ and denote by fl D the set of D n n_tuples 
with ij taking values in {1, 2, ■ • • , D}. 

The atoms of the partitions £[o,n-i] describe segments of trajectories up to time 
n encoded by the atoms of £ that are traversed at successive times. The richness 
in diverse trajectories, that is the degree of irregularity of the motion (as seen 
with the accuracy of the given coarse-graining), can be measured by the Shannon 
entropy [10J 

On the long run, £ attributes to the dynamics an entropy per unit time-step 

h (T a ,£):= lim -£^(£[0^-1]) . (17) 

n^oo Tl 

This limit is well defined p] and the Kolmogorov entropy h^{T a ) of (A*, coV, a ) 
is defined as the supremum over all finite measurable partitions [^ITfl]: 

hfj,(T a ) := sup h^,(T a , £) ■ (18) 

£ 



3.1. Symbolic Models as Classical Spin Chains 

Finite partitions £ of X provide symbolic models for the dynamical systems (X, /1, T a ) 
of Section^ whereby the trajectories {T^x}-^ are encoded into sequences {ij}j eZ 
of indices relative to the atoms E ij visited at successive times j; the dynamics cor- 
responds to the right-shift along the symbolic sequences. The encoding can be 
modelled as the shift along a classical spin chain endowed with a shift-invariant 
state [12j. This will help to understand the quantum dynamical entropy which will 
be introduced in the next Section. 

Let D be the number of atoms of a partition £ of X , we shall denote by A D 
the diagonal D x D matrix algebra generated by the characteristic functions of 
the atoms Eg and by A^' n ^ the n-fold tensor product of n copies of (Ad), that 
is the D n x D n diagonal matrix algebra A^' n_1 ' := (A D ) <S> (A D ) 1 ■ ■ ■ <g> (Ad)„_i. 
Its typical elements are of the form do ® a± • • • ® a n -\ each a-, being a diagonal 
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D x D matrix. Every A^' 9 ' := ® 9 j =p (A D )j can be embedded into the infinite tensor 
product A%:= &™ =0 (A D ) k as 

(1) <g> • • ■ ® (l) p _! ® (A D ) P ® • • ■ ® (A D ), ® (l) ?+ i ® (l) g+2 (81 • • • (19) 

The algebra A 1 ^ is a classical spin chain with a classical .D-spin at each site. 

By means of the discrete probability measure {pi}^^, one can define a com- 
patible family of states on the "local" algebras A^'™ -1 ^: 



pf n 11 (ao ® • ■ ■ ® Ctn-l) = ^ Pi ( a o)i io • ■ ■ (On-l)tn-lin-l • ( 20 ) 



Indeed, let pfiV denote the restriction to a subalgebra N C M" of a state p on 
a larger algebra M. Since ^ _ Pi ii---i n -i = / i iou---in-2J w hen n varies the local 
states pg' n X ^ are such that p^ ,n 1 [ A^' n 2 ^ = p^' n 2 ' and define a "global" state 
P£ on A^ such that pg [ A^' n ^ = p^'™ 1 . 

From the T a -invariance of p it follows that, under the right-shift o : A^ \— > A^ 1 , 

a (A^ ] ) = A% +1 ' q+1] , (21) 
the state ps of the classical spin chain is translation invariant: 

PS O CT (ao ® • • " ® O n -l) = P£ ((l)o ® (Oo)l ® ■ ■ • ® (On-l)n) ^ 

= PS {do <8> • • • ® On-l) • 

Finally, denoting by |j) the basis vectors of the representation where the matrices 
a G Ad are diagonal, local states amount to diagonal density matrices 

p pn-i] = p. \i Q )(i \ | <«! | <g> • • • <g) |z n _i)(z n _i| , (23) 

and the Shannon entropy (|TK|) to the Von Neumann entropy 

Wo.n-1]) = -Tr U°' n - 1] logp' '^ 1] l =: H, [S^] . (24) 



4. ALF-Entropy 

From an algebraic point of view, the difference between a triplet (M.,u>, G) describ- 
ing a quantum dynamical system and a triplet (Ax, c^, G Q ) as in Definition 12.11 is 
that uj and G are now a G-invariant state, respectively an automorphism over a 
non-commutative (C* or Von Neumann) algebra of operators. 
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Remark 4.1 

In finite dimension D, M. is the full matrix algebra of D x D matrices, 
the states u are given by density matrices p w , such that u(X) := Tr(p^X), 
while the reversible dynamics is unitarily implemented: Q{X) = 
UXU*. 

The quantum dynamical entropy proposed in [12] by Alicki and Fannes, ALF 
entropy for short, is based on the idea that, in analogy with what one does for 
the metric entropy, one can model symbolically the evolution of quantum systems 
by means of the right shift along a spin chain. In the quantum case the finite- 
dimensional matrix algebras at the various sites are not diagonal, but, typically, 
full matrix algebras, that is the spin at each site is a quantum spin. 

This is done by means of so called partitions of unit, that is by finite sets 
y = \ yi,y 2 ,...,y D } of operators in a O-invariant subalgebra A4o E M. such that 

D 



5>^ =1 > ( 25 ) 



1=1 



where y* denotes the adjoint of j/-. With y and the state u one constructs the 
D x D matrix with entries ^(y*^ ); such a matrix is a density matrix p\y\. 

p[y] itj := u{y* Vi ) ■ (26) 

It is thus possible to define the entropy of a partition of unit as (compare (J23J): 

^[^:=-Tr(p[y]logp[^]) • (27) 

Further, given two partitions of unit y = (yo, y±, . . . , yo \ , Z = (zo, z±, . . . , zb^J , 
of size Z), respectively B, one gets a finer partition of unit of size BD as the set 

y o Z:= (y z , ■■■ , y z B ; yiz , ■■■ , yiZ B ; ■ • • ; z/d^o, • • • , VdZb) ■ (28) 

After j time-steps, y evolves into &(y) := |@-' (yi), Q^^), • ■ ■ , (?/£>) |. Since 

9 is an automorphism, Qi (y) is a partition of unit; then, one refines Qi(y), 
^ j ^ n — 1, into a larger partition of unit 

y[o.»-i] := e«-i(y) o e n ~ 2 (y) o ... o Q(y) o y. (29) 

We shall denote the typical element of ^I -" -1 !] by 

[y l0 ' n -\ = (y in J © n ' 2 (y in J ■ ■ ■ I/*- (so) 
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Each refinement is in turn associated with a density matrix py' n ^ := p [3^ 0,n 
which is a state on the algebra M^' n_1 ' : 



»n-l, 



y 



[0,71-1] 



e=0 (M.£ ) ) l , with entries 



Moreover each refinement has an entropy 



y 



[0,n-l] 



-Trip 



y 



[0,n-l] 



logp 



The states p^' n ^ 



-y[0,n-l] 
[0,n-2] 



(31) 



(32) 

and define a global 



are compatible: p y ' n \ M. [ D ' n = p y 
state py on the quantum spin chain M^ 1 := Cg>£L (M£))£. 

Then, as in the previous Section, it is possible to associate with the quantum 
dynamical system 9) a symbolic dynamics which amounts to the right- 

shift, a : (M fl ) f i— y (M D ) £+1 , along the quantum spin half-chain (compare (|2~T| ). 
Non-commutativity makes py not shift-invariant, in general ^2j. In this case, 
the existence of a limit as in (fTTj) is not guaranteed and one has to define the 
ALF-entropy of (M,uj,Q) as 



K%M--= ^p K% (e,y) 

ycMo 



where 



h^ L M (<d,y)--= limsup-tf, 



y 



[0,n-l] 



(33a) 
(33b) 



Like the metric entropy of a partition £, also the ALF-entropy of a partition of 
unit y can be physically interpreted as an asymptotic entropy production relative 
to a specific coarse-graining. 

Remark 4.2 



The ALF-entropy reduces to the Kolmogorov metric entropy on clas- 
sical systems. This is best seen by using an algebraic characteri- 
zation of (X,p,T a ) by means of the Von Neumann algebra M.x = 
(X) of essentially bounded functions on X [201 • The character- 
istic functions of measurable subsets of X constitute a *subalgebra 
M.q C M.x\ moreover, given a partition 8 of X, the characteristic func- 
tions eE e of its atoms E#, Zg = {e^, ■ ■ • ,e£ D } is a partition of unit 
in A4 . From Q it follows that Q J a {eE e ) = e T~ j {E t ) an( ^ fr° m that 

p[zf' n ~ 1] ] = Sijpi (see JH3)), whence H^zf^^] = (S [o,n-i]) 

J hj 

(see (fT6|) and (|27|) V In such a case, the limsup in (|33bj) is actually 

a true limit and yields (fTTj) . In [TJ], the same result is obtained by 

means of the algebra Ax and of the *subalgebra W exp of exponential 

functions. 
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4.1. ALF-Entropy for Discretized (X,n,T a ) 

We now return to the classical systems (Ax, w^, Q a ) of Section [2j For later use, 
we introduce the following map defined on the torus T 2 ([0,iV) 2 ), namely [0,N) 2 
(mod N), and on its subset (Z/7VZ) 2 : 

T 2 ([0,iV) 2 ) 3x^U a (x):=NT a (^j GT 2 ([0,iV) 2 ) (34) 

The use of the *morphisms Jn,oo and Jca,N introduced in Section 12.11 makes it 
convenient to define the discretized versions of (Ax, Up, @q) as follows: 

Definition 4.1 

A discretization of (Ax, w^, Q a ) is the triplet (Dn? (C) , con 2 , ©a) where: 

D N 2 (C) is the abelian algebra of diagonal matrices acting on C^ 2 . 
u>n 2 is the tracial state given by the expectation: 

D N 2 (C) 9 M h-> uj N 2 (M) := Tr (M) . (35) 

a is the *automorphism of D N 2 (C) defined by: 

D N 2 (C) 3 M i > Q a (M) := M u a W,u aW \t) (t\ ■ (36) 

e.e(%/NZ) 2 

Remarks 4.3 

i. The expectation ui N 2 (Jn,oo (/)) corresponds to the numerical calcu- 
lation of the integral of / realized on a N x N grid on T 2 . 

ii. Q a is a *automorphism because the map (Z/7VZ) 2 3 1 i — >■ {7q, (£) is 
a bijection. For the same reason the state u^ 2 is ©^-invariant. 

iii. One can check that, given / G Ax, 

& a (JN,oo(f)):= £ /(^) I *><*|- (37) 

iv. Also, ©£ o Jiv i00 = J^ i00 o e£ for all j € Z. 

v. On the contrary, for j G IN, 9^ o ^,00 7^ <7n,oo ©i for the Sawtooth 
Maps, that is when a G" Z. 
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The automorphism Q a can be rewritten in the more familiar form 

Q a (X)= J2 x u a (t),u ai e)\£)(£\ 



c/~ 1 (s)e(z/Arz) 2 

(see Remark |Ui and ii) = U aN ^ X s>s \s) (s\\ U* >N (38) 

\ all equiv. / 
> classes ' 

= U a>N X U*, N , (39) 
where the operators U a ^ era defined by 

n N 2 3\i)^U^ N \t):=\U- 1 {i)) ■ (40) 

Remarks 4.4 

i) All of T a , T~\ T* and (T" 1 )* belong to SL 2 (Z/NZ); in particular 
these matrices are automorphisms on (Z/NZ) 2 so that, in (|38l) . one 
can sum over the equivalence classes. 

ii) The same argument as before proves that the operators in (|4*fl|) are 
unitary which is equivalent to saying that Q a is a *automorphism. 

In order to construct the ALF-entropy, we now seek a useful partition of unit in 
(D N 2,lun2, G a ); we do that by means of the subalgebra Woo ^ A* in equations 



y:=\y 3 \ = {-?= exp (27rzr,- • x) \ , (41) 



3=1 j j=1 

where j^j =: Ac (Z/NZ) 2 - (42) 
Definition 4.2 

Given a subset A of the lattice consisting of the points {rv, } as in P^Jl . 
we shall denote by y the partition of unit in (D N 2,uj N 2, Q a ) given by: 

D ( ^D ( y ^D 



j=l k ) 3=1 K V D ) j=1 



y={^)_--= (^.oo (W) !> \ -7= U'irj)} . (43) 
with W(rj) defined in 
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From the above definition, the elements of the refined partitions in (|30l) take the 
form: 



y 



[0,71-1] 



1 1 

ND% 



E 



e N 



[n n _ 1 -U2- 1 (£)+-+r il -U a (£)+T iQ -£] i_gv ^| 



(44) 



Then, the multitime correlation matrix p^~' n ^ in Q31j) has entries: 



3> 



[0,n-l] 



. . iV 2 L> n ^ 

% <3 £e(Z/NZ) 2 



■liti 

N 



E (T-i p -T-i p )-J7S(^) 



p=0 



(*\9t(n)) (g t {n)\j) 

£€(Z/NW,) 2 



K (£) = t (45) 



(46) 



with {i\g e (n)): 



1 J_ JV 

[0,71-1 



n-l 

E 

e 



e C 1 



(47) 



The density matrix p~' can now be used to numerically compute the ALF 
entropy as in (|33|) : however, the large dimension (D n x D n ) makes the computa- 
tional problem very hard, a part for small numbers of iterations. Our goal is to 
prove that another matrix (of fixed dimension N 2 x N 2 ) can be used instead of 

[0,71-1] 

Py ■ 

Proposition 4.1 

Let Q (n) be the iV 2 x N 2 matrix with entries 

& llfe (n):= (^(n)|^(n)> (48) 

given by the scalar products of the vectors \g^{n)) G Hd™ = C 15 ™ 
in (JITJl . Then, the entropy of the partition of unit 3^°' n_1 ' with ele- 
ments (jlHJ) is given by: 



V2 



3> 



[0,71-1] 



Proof: £ (n) is hermitian and from ijlTJ) it follows that Ty Hn2 Q 



n 



(49) 



I. 



Let Hn 2 = C , H:= Hd™ ® H-n 2 and consider the projection /tty = IVO (VH onto 



£e(z/ivz) 2 



(n))<8)|£). 



(50) 



We denote by S x the restriction of p^ to the full matrix algebra M\\= M D n (C) 
and by S 2 the restriction to M 2 := (C). It follows that: 

Tr HDn (Ei • mi) = (^|mi ® 1 2 = (^1 m i l&) > Vm i e M i ' 

•te(Z/7VZ) 2 
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Thus, from ijlBjl . 



Si = P^ 1] = £ I I • (51) 

£e(z/7vz) 2 



On the other hand, from 

Tr Wjv2 (£ 2 ■ m 2 ) = {ip\ li® m 2 

Yl ^12 ( n )\9£i (n)) (4|m 2 |£i) , Vm 2 G M 2 , 
e 1 ,i 2 e(i./Nz) 2 

it turns out that £ 2 = (n), whence the result follows from Araki-Lieb's inequal- 
ity |2H ■ 

We now return to the explicit computation of the density matrix Q (n) in Propo- 
sition HUJ By using the transposed matrix T^ r , the vectors ijlTJ) now read 

(i|^(n)> = ^ f e-^ir« (52) 

n-l 

(mod AT) (53) 

p=0 

where we made explicit the various dependencies of (fBTfl) on n the time-step, N 
the inverse lattice-spacing, the chosen set A of t\,-'s and the a parameter of the 
dynamics in SL 2 (Z/NZ) 2 . 

In the following we shall use the equivalence classes 

[r] := [i G | f^f (<) = r e (Z/iVZ) 2 (mod N)} , (54) 

their cardinalities # [r] and, in particular, the frequency function v^lf 1 

(Z/NZ) 2 3 r _ vflf (r) := . (55) 

Proposition 4.2 

The Von Neumann entropy of the refined (exponential) partition of 
unit up to time n — 1 is given by: 



TT 



y 



[0,n-l] 



J2 ^ N (r)\ogv^ N (r) (56) 

r£(Z/NW,) 2 
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Proof: Using (f52l) . the matrix Q (n) in Proposition 14. II can be written as: 

0(») = ^E <£( n ) I ' ^ 

{Z\k{n)) = -e— f ^ (58) 
The vectors | /< (n) ) G 7^ = C^ 2 are such that (f t (n)| (n)) = 5 { V )iN (n) „ 

J A, a W ' J A, a UJ 

where with 5^ is the A^-periodic Kronecker delta. For sake of simplicity, we say 
that | fi (n) ) belongs to the equivalence class [r] in (JKIJ) if * G [r] ; vectors in differ- 
ent equivalence classes are thus orthogonal, whereas those in a same equivalence 
class [r] are such that 



\i€[r] ) ie[r] 



D n ^ N (r) (i^eir)) (e(r)\£ 2 ) 



(£\e(r)) = ^^EH N2 
Therefore, the result follows from the spectral decomposition 
S(A= £ ^ N (r)\e(r))(e(r)\ 

re(W,/NZ) 2 



5. Analysis of entropy production 

In agreement with the intuition that finitely many states cannot sustain any lasting 
entropy production, the ALF-entropy is indeed zero for such systems jT2| - How- 
ever, this does not mean that the dynamics may not be able to show a significant 
entropy rate over finite interval of times, these being typical of the underlying 
dynamics. 

As already observed in the Introduction, in quantum chaos one deals with 
quantized classically chaotic systems; there, one finds that classical and quantum 
mechanics are both correct descriptions over times scaling with \ogh~~ 1 . Therefore, 
the classical-quantum correspondence occurs over times much smaller than the 
Heisenberg recursion time that typically scales as h~ a , a > 0. In other words, for 
quantized classically chaotic systems, the classical description has to be replaced 
by the quantum one much sooner than for integrable systems. 

In this paper, we are considering not the quantization of classical systems, 
but their discretization; nevertheless, we have seen that, under certain respects, 
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quantization and discretization are like procedures with the inverse of the number 
of states N playing the role of h in the latter case. 

We are then interested to study how the classical continuous behaviour emerges 
from the discretized one when N — > oo; in particular, we want to investigate the 
presence of characteristic time scales and of "breaking-times" tb, namely those 
times beyond which the discretized systems cease to produce entropy because their 
granularity takes over and the dynamics reveals in full its regularity. 

Propositions 14.11 and 14.21 afford useful means to attack such a problem nu- 
merically. In the following we shall be concerned with the time behavior of the 
entropy of partition of units as in Definition 14.21 the presence of breaking-times 
tb (A,N,a), and their dependence on the set A, on the number of states N and 
on the dynamical parameter a. 

As we shall see, in many cases tb depends quite heavily on the chosen partition 
of unit; we shall then try to cook up a strategy to find a tb as stable as possible upon 
variation of partitions, being led by the idea that the "true" tb has to be strongly 
related to the Lyapounov exponent of the underlying continuous dynamical system. 

Equations (|4*9|) and (|56|) allow us to compute the Von Neumann entropy of 
the state p~' n , if we were to compute the ALF-entropy according to the defini- 
tions (|33|) . the result would be zero, in agreement with fact that the Lyapounov 
exponent for a system with a finite number of states vanishes. Indeed, it is suf- 



ficient to notice that the entropy H u 



• N 2 



y[0,n- 



is bounded from above by the 



entropy of the tracial state jt^^n 2 - 



that is by 21ogA^; therefore the expression 



hw N2 ,w x {oi, A,n) 



n 



N 2 



y 



[0,71.-1] 



(59) 



goes to zero with n — > 0. It is for this reason that, in the following, we will focus 
upon the temporal evolution of the function h^^^^a, A, n) instead of taking its 
limsup over the number of iterations n. 

In the same spirit, we will not take the supremum of (|59|) over all possible 
partitions y (originated by different A); instead, we will study the dependence 



of h„ 



N 2 



(a, A,n) on different choices of partitions. In fact, if we vary over all 



possible choices of partitions of unit, we could choose A = (Z/ATZ) in (|4*2|) . that 
is D = N 2 ; then summation over all possible r G (Z/A^Z) 2 would make the matrix 



elements {n) in (|4*8|) equal to 



N 2 



whence H„ 



N 2 



y[0,n-l] 



2 log 
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5.1. The case of T a e GL 2 (T 2 ) 

The maximum of 2 is reached when the frequencies 

v^f : (Z/iVZ) 2 - [0, 1] 

become equal to 1/N 2 over the torus: we will see that this is indeed what happens 
to the frequencies with n — ► oo. The latter behaviour can be reached in 

various ways depending on: 

• hyperbolic or elliptic regimes, namely on the dynamical parameter a; 

• number of elements (D) in the partition A; 

• mutual location of the D elements T{ in A. 

For later use we introduce the set of grid points with non-zero frequencies 

£e(Z/NZ) 2 , ^ N (£)^o) ■ (60) 



r (n),N. 
1 A,o •" 



N 



5.1.1. Hyperbolic regime with D randomly chosen points r» in A 

In the hyperbolic regime corresponding to a G Z \ {—4,-3,-2,-1, 0}, 
tends to increase its cardinality with the number of time-steps n. Roughly speak- 
ing, there appear to be two distinct temporal patterns: a first one, during which 

# (r^ a N J ~ D n ^ A^ 2 and almost every ~ D~ n , followed by a second one 

characterized by frequencies frozen to v£]f (£) = ^y, W G (Z/A^Z) 2 . The sec- 
ond temporal pattern is reached when, during the first one, T^ a N has covered the 
whole lattice and D n ~ A^ 2 . 

From the point of view of the entropies, the first temporal regime is charac- 
terized by 

H u, N 2 (a, A, n) ~ n ■ log D , K N2iWao (a, A,n) ~ log£> , 
while the second one by 

2 log AT 

(a, A, n) ~ 2 log AT , K N2>Woo (a, A, n) — ■ 

The transition between these two regimes occurs at n — \og D N 2 . However this 
time cannot be considered a realistic breaking-time, as it too strongly depends on 
the chosen partition. 
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Figure [2] (columns a and c) shows the mechanism clearly in a density plot: 
white or light-grey points correspond to points of T\^ N and their number increases 
for small numbers of iterations until the plot assume a uniform grey color for large n 

The linear and stationary behaviors of H w 2 (a, A, n) are apparent in fig. El 
where four different plateaus (2 log N) are reached for four different N, and in fig. El 
in which four different slopes are showed for four different number of elements in 
the partition. With the same parameters as in fig. El fig.EJshows the corresponding 
entropy production h u JV 2,w 00 ( Q! ) A, n). 

5.1.2. Hyperbolic regime with D nearest neighbors T{ in A 



In the following, we will consider a set of points A = {ri} i=1 D very close to each 
other, instances of which are as below: 




D = 5 D = A D = 3 D = 2 



Figure 1: Several combinations of D nearest neighbors in A for different values D. 



From eqs. (|54H55l) . the frequencies (£) result proportional to how many 

strings have equal images £, through the function /a a in (|53|) . Due to the fact 

that [T a } u = [T a ] 21 = 1, non-injectivity of /^a occurs ver Y frequently when 
{ri} are very close to each other. This is a dynamical effect that, in continuous 
systems [H], leads to an entropy production approaching the Lyapounov exponent. 
Even in the discrete case, during a finite time interval though, h WN2 ^ x (a, A,n) 
exhibits the same behavior until H w a reaches the upper bound 21ogA^. From 
then on, the system behaves as described in subsection I5.1.1| and the entropy 
production goes to zero as: 

K n2 ,Woo (a, A, n) ~ i (see fig. EJ) . 

Concerning figure El (column d), whose corresponding graph for h^^^^^l, A,n) 
is labeled by > in fig. we make the following consideration: 

• for n = 1 the white spot corresponds to five r,i grouped as in fig. [TJ 
In this case K 2tWoo (l, A, 1) = logD = log 5; 
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• for n G [2, 5] the white spot begins to stretch along the stretching direction 
of T\. In this case, the frequencies i/^ are not constant on the light-grey 
region: this leads to a decrease of h w ^^^(1, A, n); 

• for n G [6, 10] the light-grey region becomes so elongated that it starts feeling 
the folding condition so that, with increasing time-steps, it eventually fully 
covers the originally dark space. In this case, the behavior of h u zyy^l, A, n) 
remains the same as before up to n = 10; 

• for n — 11, T^ a N coincides with the whole lattice; 

• for larger times, the frequencies v^\ N tend to the constant value on 
almost every point of the grid. In this case, the behaviour of the entropy 
production undergoes a critical change (the crossover occurring at n — 11) 
as showed in fig. 

Again, we cannot conclude that n — 11 is a realistic breaking-time, because once 
more we have strong dependence on the chosen partition (namely from the number 
D of its elements). For instance, in fig.0 one can see that partitions with 3 points 
reach their corresponding "breaking-times" faster than that with D = 5; also they 
do it in an A^-dependent way. 

For a chosen set A consisting of D elements very close to each other and A" very 
large, h w 2 yv^a. A, n) ~ log A (which is the asymptote in the continuous case) 
from a certain n up to a time tb- Since this latter is now partition independent, 
it can properly be considered as the breaking-time of the system; it is given by 

r B = \og x N 2 - (61) 



It is evident from equation (|fiTl) that if one knows tb 
then also log A is known. Usually, one is interested 
in the latter which is a sign of the instability of the 
continuous classical system. In the following we develop 
an algorithm which allows us to extract log A from 
studying the corresponding discretized classical system 
and its ALF-entropy. 



In working conditions, A^ is not large enough to allow for n being smaller than tb\ 
what happens in such a case is that Hu n2 (, a > A, ra) ~ 2 log A^ before the asymptote 
for h^^^^^a, A, n) is reached. Given /i Wjv2iWoo (a, A, n) for n < t b , it is thus 
necessary to seek means how to estimate the long time behaviour that one would 
have if the system were continuous. 
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Remarks 5.1 

When estimating Lyapounov exponents from discretized hyperbolic 
classical systems, by using partitions consisting of nearest neighbors, 
we have to take into account some facts: 



a. h^^^^^a, A,n) does not increase with n; therefore, if-D < A, ^^.Woo 
cannot reach the Lyapounov exponent. Denoted by logA(-D) the 
asymptote that we extrapolate from the data a , in general we have 
A (D) ^\ogD < A. For instance, for a = 1, A = 2.618 ... > 2 and 
partitions with D = 2 cannot produce an entropy greater then log 2; 
this is the case for the entropies below the dotted line in figs. EleEJ 

b. partitions with D small but greater than A allow log A to be reached 
in a very short time and A (D) is very close to A in this case; 

c. partitions with D ^> A require very long time to converge to log A 
(and so very large N) and, moreover, it is not a trivial task to deal 
with them from a computational point of view. On the contrary the 
entropy behaviour for such partitions offers very good estimates of A 
(compare, in fig. > with o, a, o and □) ; 

d. in order to compute A (and then tb, by (|HT|) ). one can calculate A (D) 
for increasing D, until it converges to a stable value A; 

e. due to number theoretical reasons, the UMG on (Z/NZy present 
several anomalies. An instance of them is showed in fig. El (col. f), 
where a partition with five nearest neighbors on a lattice of 200 x 200 
points confines the image of f^ N (under the action of a T a map 
with a = 17) on a subgrid of the torus. In this and analogous cases, 
there occurs an anomalous depletion of the entropy production and no 
significant information is obtainable from it. To avoid this difficulties, 
in Section E21 we will go beyond the UMG subclass considered so far. 



5.1.3. Elliptic regime (a E {-1, -2, -3}) 

One can show that all evolution matrices T a are characterized by the following 
property: 

Tl = aT a -\ , a:={a + 2) • (62) 
a ^w JV 2.w 00 ( a i A, n) may even equal logA(D) from the start. 
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In the elliptic regime a G {—1, —2, —3}, therefore a G {—1, 0, 1} and relation (|62l) 
determines a periodic evolution with periods: 

T_l = -t (Tj = 1) (63a) 

T_ 2 2 = -1 (T_ 2 4 = 1) (63b) 

Tj = +1 • (63c) 

It has to be stressed that, in the elliptic regime, the relations (ff)3|) do not hold 
"modulo N", instead they are completely independent from N. 

Due to the high degree of symmetry in relations (jfilMvfl) . the frequencies i^j^ N 
are different from zero only on a small subset of the whole lattice. 

This behavior is apparent in fig. [2]: col. b, in which we consider five randomly 
distributed 7*j in A, and in fig. El: col. e, in which the five are grouped as in fig.EJ 
In both cases, the Von Neumann entropy H u 3 (n) is not linearly increasing with 
n (see figure [HJ), instead it assumes a log-shaped profile ( up to the breaking-time, 
see fig. BJ. 

Remark 5.2 

The last observation indicates how the entropy production analysis 
can be used to recognize whether a dynamical systems is hyperbolic or 
not. If we use randomly distributed points as a partition, we observe 
that hyperbolic systems show constant entropy production (up to the 
breaking-time), whereas the others do not. 

Moreover, unlike hyperbolic ones, elliptic systems do not change their 
behaviour with N (for reasonably large N) as clearly showed in fig. 0J 
in which elliptic systems (a = —2) with four different values of N give 
the same plot. On the contrary, we have dependence on how rich is the 
chosen partition, similarly to what we have for hyperbolic systems, as 
showed in fig. 



5.1.4. Parabolic regime (a G {0,4}) 

This regime is characterized by A = A -1 = ±1, that is log | A | = (see Remark I2~T| 
a). These systems behave as the hyperbolic ones (see subsections 15. 1.11 and I5.1.2B 
and this is true also for the the general behavior of the entropy production, apart 
from the fact that we never fall in the condition (a.) of Remark 15.11 Then, 
for sufficiently large N, every partition consisting of D grouped will reach the 
asymptote log | A | =0. 
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5.2. The case of Sawtooth Maps 

The Sawtooth Maps [TH|fTH]. are triples (X,fj,,S a ) where 

X = T 2 = R 2 /Z 2 = {x = (x 1 ,x 2 ) (mod 1)} (64a) 

d^(cc) = da;i da; 2 , (64c) 

where (•) denotes the fractional part of a real number. Without (■), f)64b[) is not 
well defined on T 2 for not-integer a; in fact, without taking the fractional part, the 
same point x = x + nET 2 ,nEl> 2 , would have (in general) S a (x) ^ S a (x + n) . 
Of course, (•) is not necessary when a G Z. 

The Lebesgue measure defined in (|64cj) is invariant for all a G E. 

After identifying x with canonical coordinates (q,p) and imposing the (mod 1) 

condition on both of them, the above dynamics can be rewritten as: 

q ', =q+ / 1 \ ( m0d1 )' ( 65 ) 
p =p + a(q) 

This is nothing but the Chirikov Standard Map [I] in which — ^ sin(27rg) is re- 
placed by (q). The dynamics in (|65|) can also be thought of as generated by the 
(singular) Hamiltonian 

H(q,p,t) = ^-a^f5 p (t), (66) 

where 5 p (t) is the periodic Dirac delta which makes the potential act through 
periodic kicks with period 1. 

Sawtooth Maps are invertible and the inverse is given by the expression 

^(3-U0(C"0(3> (modl) (67) 

or, in other words, 

(mod 1) . (68) 

p = —a q + p 

It can indeed be checked that S a (S^ 1 (x)) = S^ 1 (S a (x)) = (x) , Va? G T 2 . 



Remarks 5.3 
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i. Sawtooth Maps {S a } are discontinuous on the subset 

7 := {x = (0,p) , pel} G T 2 : two points close to this border, 
A:= (e,p) and B:= (1 — £,p), have images that differ, in the e — > 
limit, by a vector c?^ (A, B) = (a, a) (mod 1). 

ii. Inverse Sawtooth Maps {S^ 1 } are discontinuous on the subset 
7_i := 5* Q (70) — {x — (p,p) , p G T} G T 2 : two points close to this 
border, A := (p + e, p — e) and B:= (p — e, p + e) , have images that 
differ, in the e — > limit, by a vector d^l 1 (A, B) = (0, a) (mod 1). 

'-'a 

iii. The hyperbolic, elliptic or parabolic behavior of Sawtooth maps is 
related to the eigenvalues of ( 1 + a \ ) exactly as in Remark EHJii. 

iv. The Lebesgue measure in (|64cjl is S^-invariant. 

From a computational point of view, the study of the entropy production in the 
case of Sawtooth Maps S a is more complicated than for the T a 's. The reason to 
study numerically these dynamical systems is twofold: 

• to avoid the difficulties described in Remark [5. II (e.): 

• to deal, in a way compatible with numerical computation limits, with the 
largest possible spectrum of accessible Lyapounov exponent. We know that 
for « 6 Zf| {non elliptic domain}, 

In order to fit log X a ( log X a being the Lyapounov exponent corresponding 
to a given a) via entropy production analysis, we need D elements in the 
partition (see points b. and c. of Remark 15.11) with D ^ A a . Moreover, 
if we were to study the power of our method for different integer values of 
a we would be forced forced to use very large D, in which case we would 
need very long computing times in order to evaluate numerically the entropy 
production h u 2 ^(q, A, n) in a reasonable interval of times n. Instead, 
for Sawtooth Maps, we can fix the parameters (N, D, A) and study A Q for 
a confined in a small domain, but free to assume every real value in that 
domain. 



In the following, we investigate the case of a in the hyperbolic regime with D 
nearest neighbors in A, as done in subsection 15. 1.21 In particular, figures |9HT2l) 
refer to the following fixed parameters: 
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N = 38 ; n max =5 ; D = 5 



A 




a : from 0.00 to 1.00 with an incremental step of 0.05. 



First, we compute the Von Neumann entropy (|49|) using the (hermitian) matrix 
Gii& ( n ) defined in ijlHjl . This is actually a diagonalization problem: once that 
the N 2 eigenvalues {r]i}^ =l are found, then 



N 



^rjilogrji ■ (69) 



Then, from (|59j) . we can determine ^a; 2,W 00 ( Q! >A, n). In the numerical example, 
the (A-dependent) breaking-time occurs after n = 5; for this reason we have 
chosen n max = 5. In fact, we are interested in the region where the discrete system 
behaves almost as a continuous one. 

In figurelHl the entropy production is plotted for the chosen set of a's: for very 
large N (that is close to the continuum limit, in which no breaking-time occurs) 
all curves (characterized by different a's) would tend to logA a with n. 

One way to determine the asymptote logA^ is to fit the decreasing function 
h^ N2t Wa ( a :-^J n ) over the range of data and extrapolate the fit for n — > oo. Of 
course we can not perform the fit with polynomials, because every polynomial 
diverges in the n — > oo limit. 

A better strategy is to compactify the time evolution by means of a isomorphic 
positive function s with bounded range, for instance: 

2 

IN* 9 n i — > s n := - arctan (n - 1) G [0, 1] ■ (70) 

7T 

Then, for fixed a, in fig. HH] we consider n max points (s n , h w ~ ^(a, A, n)) and 
extract the asymptotic value of ^^^(a, A,n) for n — > oo, that is the value of 

hwrfiWoo s ~ 1 (*)) f° r ^ > as follows. 

Given a graph consisting of m e {2, 3, • • • , n max } points, in our case the first 
m points of curves as in fig. [TOl namely 

{( s i , h u N2 ,Woo{®, A, 1)) , (s 2 , hu^-Wooia, A, 2)) , • • • , (s m , ^^(a, A,m))} , 
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the data are fit by a Lagrange polynomial V m (t) (of degree m — 1) 



V m it) 



(71a) 



i=i 




^ w2 ,Woc(a ; A,i) • 



(71b) 



The value assumed by this polynomial when t — > 1~ (corresponding to n — > oo) 
will be the estimate (of degree m) of the Lyapounov exponent, denoted by I™: the 
higher the value of m, the more accurate the estimate. From (I7TJ) we get: 



The various I™ are plotted in figure HU as functions of m for all considered a. The 
convergence of I™ with m is showed in figure [TJl together with the theoretical 
Lyapounov exponent logA^; as expected, we find that the latter is the asymptote 
of {/™} m with respect to the polynomial degree m. 

The dotted line in fig. HH extrapolates 21 a-curves in compactified time up to 
t — 1 using five points in the Lagrange polynomial approximation. 

6. Conclusions 

In this paper, we have considered discretized hyperbolic classical systems on the 
torus by forcing them on a squared lattice with spacing jr. We showed how the 
discretization procedure is similar to quantization; in particular, following the anal- 
ogous case of the classical limit h i — > 0, we have set up the theoretical framework 
to discuss the continuous limit iV i — ► oo. Furthermore, using the similarities be- 
tween discretized and quantized classical systems, we have applied the quantum 
dynamical entropy of Alicki and Fannes to study the footprints of classical (con- 
tinuous) chaos as it is expected to reveal itself, namely through the presence of 
characteristic time scales and corresponding breaking-times. Indeed, exactly as in 
quantum chaos, a discretized hyperbolic system can mimick its continuous part- 
ner only up to times which scale as logiV. We have also extended the numerical 
analysis from the so called Arnold cat maps to the discontinuous Sawtooth maps, 
whose interpretation in the theoretical frame set up in this work will be discussed 
in a forthcoming paper. 
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Figure 2: Density plots showing the frequencies ^ ,JV in two hyperbolic regimes 
(columns a and c) and an elliptic one (col. b), for five randomly distributed r*j in 
A with N = 200. Black corresponds to Vjfy N = 0. In the hyperbolic cases, v^f 1 
tends to equidistribute on (Z/iVZ) 2 with increasing n and becomes constant when 
the breaking-time is reached. 
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Figure 3: Density plots showing Vj^ N in two hyperbolic (columns d and f) and 
one elliptic (col. e) regime, for five nearest neighboring r« in A (N = 200). Black 
corresponds to Vj^ N = 0. When the system is chaotic, the frequencies tend to 
equidistribute on (Z/iVZ) 2 with increasing n and to approach, when the breaking- 
time is reached, the constant value Col. (f) shows how the dynamics can be 
confined on a sublattice by a particular combination (a, N, A) with a corresponding 
entropy decrease. 
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Number of iterations n 



Figure 4: Von Neumann entropy H u 2 (n) in four hyperbolic (a — 1 for o, A, o, 
□) and four elliptic (a = —2 for >) cases, for three randomly distributed r*j in A. 
Values for N are: o = 500, a = 400, o = 300 and □ = 200, whereas the curve 
labeled by > represents four elliptic systems with N e {200,300,400,500}. 




o J , , , , , , , , , , , , , V 

1 2 3 4 5 6 7 8 9 10 1 1 12 1 3 14 15 



Number of iterations n 

Figure 5: Von Neumann entropy 2 (n) in four hyperbolic (a = 1) cases, for D 
randomly distributed r*j in A, with N = 200. Value for D are: o = 5, a = 4, o = 3 
and □ = 2. The dotted line represents H w 2 (n) = log A • n where log A = 0.962 . . . 
is the Lyapounov exponent at a = 1. 
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1 2 3 4 5 6 7 8 9 10 1 1 1 2 1 3 14 1 5 
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Figure 6: Entropy production /i Wjv2 ^(a, A, n) in four hyperbolic (a = 1) cases, 
for D randomly distributed in A, with N = 200. Values for D are: o = 5, 
a = 4, o = 3 and □ = 2. The dotted line corresponds to the Lyapounov exponent 
log A = 0.962. . . at a = 1. 
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Figure 7: Entropy production /i Wjv2 ^(a, A, n) in five hyperbolic (a = 1) cases, 
for D nearest neighboring points in A. Values for (N, D) are: > = (200, 5), o = 
(500,3), a = (400,3), o = (300,3) and □ = (200,3). The dotted line corresponds 
to the Lyapounov exponent log A = 0.962 ... at a = 1 and represents the natural 
asymptote for all these curves in absence of breaking-time. 
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Number of iterations n 



Figure 8: Von Neumann entropy H u 2 (n) in four elliptic (a = —2) cases, for D 
randomly distributed T*j in A, with N = 200. Value for D are: o = 5, A = 4, o = 3 
and □ = 2. 
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Figure 9: Entropy production h u ^^^(a, A, n) for 21 hyperbolic Sawtooth maps, 
relative to a for a cluster of 5 nearest neighborings points r*j in A, with N = 38. 
The parameter a decreases from a — 1.00 (corresponding to the upper curve) to 
a = 0.00 (lower curve) through 21 equispaced steps. 
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Figure 10: The solid lines correspond to (s n , h w 2t yv x ( a , A, n)), with n G 
{1,2,3,4,5}, for the values of a considered in figure El Every a-curve is con- 
tinued as a dotted line up to (1, 1^), where l b a is the Lyapounov exponent extracted 
from the curve by fitting all the five points via a Lagrange polynomial V m (t) . 




Degree of accuracy m 



Figure 11: Four estimated Lyapounov exponents I™ plotted vs. their degree of 
accuracy m for the values of a considered in figures E3 and HH 
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Hyperbolicity parameter a 



Figure 12: Plots of the four estimated of Lyapounov exponents I™ of figure [TT] 
vs. the considered values of a. The polynomial degree m is as follows: o = 2, 
a = 3, o = 4 and □ = 5. The soli d line corr esponds to the theoretical Lyapounov 
exponent log A a = log (a + 2 + y/ a {a + 4) ) — log 2. 
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